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Abstract 

Geomagnetic effects distort the zenith angle distribution of sub-GeV and few- 
GeV atmospheric neutrinos, breaking the up-down symmetry that would be 
present in the absence of neutrino oscillations and without a geomagnetic field. 
The geomagnetic effects also produce a characteristic azimuthal dependence 
of the i^-fluxes, related to the well known east-west effect, that should be 
detectable in neutrino experiments of sufficiently large mass. We discuss 
these effects quantitatively. Because the azimuthal dependence is in first 
order independent of any oscillation effect, it is a useful diagnostic tool for 
studying possible systematic effects in the search for neutrino oscillations. 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



The flux of atmospheric neutrinos at a fixed value of the energy E v depends both on 
the zenith and azimuth angles (6 Z and if). The angular dependence originates from two, or 
possibly three sources: 

1. The development of cosmic ray showers in the atmosphere. 

2. Geomagnetic effects on the primary cosmic ray flux. 

3. Neutrino oscillations. 

The development of an hadronic shower induced by a primary particle of given energy and 
mass depends only on the zenith angle. Cascades at large zenith angle develop in a relatively 
less dense part of the atmosphere, so that decay to neutrinos is enhanced at large angle. In 
fact, the calculation of the secondary beam depends only on | cos8 z \ because a line of sight 
entering the detector from below the horizon with zenith angle 9° ut > tc/2, corresponds to a 
trajectory entering the atmosphere with 6™ = n — 9° ut . Apart from small effects due to the 
different average temperature profiles of the atmosphere at different geographical locations, 
the development of a shower does not depend on the position of its impact point on the 
earth's surface. Thus, production of secondary particles in the atmosphere is symmetric 
under the reflection cos6 z <-> — cos9 z . 

For a given energy spectrum, neutrino oscillations also depend only on zenith angle; 
however, the dependence is strongly asymmetric because the pathlengths corresponding to 
the directions ± cos 9 Z are very different. For down-going particles, the neutrino pathlengths 
are in the range from ~ 10 to ~ 500 km [p]], whereas up-going neutrinos have L ~ 10 4 km. 

Geomagnetic effects modify the spectrum of primary cosmic rays up to tens of GeV in a 
way that depends on azimuth as well as zenith. Since the neutrino flux is a convolution of 
the primary spectrum with the yield of neutrinos per primary particle, neutrinos with ener- 
gies below a few GeV carry the imprint of these geomagnetic effects. The geomagnetic field 
prevents primary cosmic rays of low rigidity from reaching the atmosphere. This suppression 
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depends on both the detector location, being lowest (highest) at the magnetic poles (equa- 
tor), and on the line of sight considered. For directions from below the horizon the effect 
must obviously be calculated for the geomagnetic field at the position where the cosmic-ray 
trajectories enter the volume of the atmosphere. The nuclear component of cosmic rays is 
positively charged, and this introduces a dependence on the azimuth angle, the celebrated 
east-west effect. The neutrino flux is highest (lowest) for directions coming from the west 
(east). 

For most interactions of atmospheric neutrinos in present detectors the direction of the 
neutrino is not fully reconstructed. Typically, the direction of the charged lepton will be used 
to indicate the direction of the event. For detected leptons with momenta in the interval 
0.2 < E^ e < 1 GeV the east-west asymmetry is of order 30%, after taking into account 
the dilution of the effect due to the broad distribution of angles between the charged lepton 
and the neutrino. An effect of this size should be readily measurable by high statistics 
experiments. For events in the multi-GeV range, when the initial neutrino energy is of 
order of several GeV the east-west effect is reduced to ~ 10%. For a detailed study it is 
necessary to consider the exact geographical location of the detector, and the interval of 
neutrino energy that is detected. 

In the presence of neutrino oscillations the zenith angle distribution of the detected events 
can be significantly deformed; however, the asymmetry in azimuth remains unchanged to 
first order because the neutrino pathlength does not depend on ip. Only to second order 
would the deformation of the neutrino energy spectrum by oscillations lead to a slight 
modification of the distribution in azimuth. The azimuthal distribution depends only on 
the filtering of the primary cosmic rays through the geomagnetic field as viewed from each 
detector. Study of the azimuthal dependence of neutrino interactions can therefore be a 
valuable diagnostic tool, both to validate calculations of the neutrino flux and to check the 
quality of detector performance, for example, to demonstrate that the determination of the 
lepton directions has the expected resolution. Moreover, since the geomagnetic effects are 
the only known mechanism (besides //-oscillations) that can produce an up-down asymmetry 
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for the neutrino fluxes, a measurement of the east-west effect for neutrinos would establish 
the size of the geomagnetic effects, and greatly help in limiting the possible importance of 
geomagnetic effects on the zenith angle distributions. The existence of distortions of the 
zenith angle distributions due to neutrino oscillations could then be more clearly identified. 
This discussion can be summarized in the following equation: 

<fi Ua (E u ,tt u ,x d ) = J2<fiA(E ) F M [p /Z,tt ,x(tt )] x 

A 

Y,-T^( E ^ A,E ,\cose z \) (P u ^ Ua (E u ,cos6 z ,{m*,U aj })) (1) 

where <p Va (E v ,Q l ,,x c i) is the flux of neutrinos of flavor a with energy E v and direction Q v 
observable in a detector located at a position Xj. 4>a(E ) is the flux of primary cosmic rays of 
energy Eq in the vicinity of the earth, but at a distance r sufficiently large so that the effects 
of the geomagnetic field are negligible (r > 10 r e , where r e is the radius of the earth). This 
flux is isotropic to a very good approximation. The cutoff factor Fm takes into account the 
effects of the geomagnetic field. It depends on the rigidity R = po/Z and direction Q of the 
primary particle and on the position x in where its trajectory first intersects the atmosphere. 
dn Vp /dE u is the average number of neutrinos of flavor j3 produced by a primary particle of 
mass A and energy Eq. It depends on the mass, energy and zenith angle of the primary 
particle. Finally, (P v „-+ Va ) is the oscillation probability for the transition v@ — > u a averaged 
over the position of creation of the neutrinos. It depends on the energy and trajectory of 
the neutrino and on the oscillation parameters; that is, the mass eigenvalues rrij and the 
mixing matrix U a j that relates mass and flavor eigenstates. 

We first discuss the calculation of the probability of penetration through the geomagnetic 
field. Then in §3 we describe the Monte Carlo convolution expressed in Eq. J. Results and 
discussion follow. 
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II. GEOMAGNETIC EFFECTS 



To a first approximation the effects of the field can be described simply by a cutoff 
rigidity R c (xd, ft), which is a function of the detector position x and the direction Q, such 
that all rigidities smaller (larger) than R c are forbidden (allowed); that is: 

F M (R,n,£ d ) = Q[R-R D (9 s ,<p,£ d )], (2) 

where Q(x) is the Heavyside function (Q(x) — for x < 1, Q(x) = 1 for x > 1). 

In the case of a dipolar magnetic field that fills the entire space, it is possible to compute 
the cutoff rigidity exactly: 

R c = R s (r, A M , g , <p) = (^) ( 77—77- T\ XM a • 1 ' ( 3 ) 

\2rv [ [1 + (1 + cos Am sin u z sin ip) l i z \ z J 

where M is the magnetic dipole moment, r (the distance from the dipole center) and Am 
(the magnetic latitude) describe the detector position, 8 Z is the zenith angle and (p is an 
azimuthal angle, with if = (|) indicating the north (west) direction. For the earth 
M ~ 8.1 x 10 25 Gauss cm 3 , which corresponds to a polar magnetic field of 0.62 Gauss. 
The quantity M/(2r 2 e ) ~ 59.4 GV corresponds to the rigidity of a particle in a circular orbit 
of radius r e in the earth's magnetic equatorial plane. 

Stormer's formula (Eq. ^) gives a good idea of the magnitude of the geomagnetic cutoffs, 
but it has limited accuracy because the geomagnetic field is only approximately an offset 
dipole. The formula also generally underestimates the cutoffs because it neglects the shadow 
of the Earth, i.e. allows the penetration of charged particles with trajectories that would 
have intersected the surface of the Earth. More exact calculations can be done using the 
backtracking technique |§ and more realistic models of the geomagnetic field 0. 

In the backtracking technique, to establish if a particle with charge Z and momentum p 
traveling from interplanetary space can reach a final point x close to the surface of the earth 
arriving from the direction Q, one integrates the equation of motion for a particle with op- 
posite charge and reflected momentum starting from this final position. If the backtracked 
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anti-particle reaches infinity, we can assume that the trajectory is allowed (Fm = 1), if 
the backtracked particle is trapped in the geomagnetic field or if its trajectory intersects 
the surface (r = r®) the trajectory is considered forbidden (Fm = 0). Such a calculation 
was performed in Ref. [|J, considering as 'trapped' those trajectories that remained confined 
within r < 30 r ffi for a pathlength longer / > 500 r®. Instead of the sharp rigidity cutoff 
predicted by equation [3|, above which all particles from a particular direction reach the at- 
mosphere, one ecounters a quite different situation: in the vicinity of Rs particle trajectories 
change rapidly with the rigidity, and the sharp cutoff is replaced with a series of allowed 
(Fm = 1) and disallowed (Fm = 0) rigidities - the penumbra region. 

The function F M used in this paper is calculated by backtracking particles for a set of 
rigidities at A cos9 z of 0.02 and A0 of 5°. The results are than averaged for angular bins of 
Acos# z = 0.1 and A0 of 30°, i.e. using the cutoffs for 36 directions for every rigidity value. 
The cutoffs are thus replaced with the probability per angular bin for a cosmic ray of given 
rigidity to reach vertical altitude of 20 km and interact in the atmosphere. 

III. MONTECARLO CALCULATION 

The calculation of the neutrino flux that we use here is described in Refs. |5]|| . Yields of 
neutrinos are calculated separately for a grid of energies for primary protons and neutrons. 
In the energy range of interest here, approximately 80% of the incident nucleons are free 
protons. Most of the rest are neutrons and protons in primary alpha particles. For each 
direction (20 bins of cos 9 Z and 12 bins of <f>) the yields are folded with the primary spectrum 
to obtain the neutrino flux. The primary spectrum is weighted with the cutoffs averaged 
over the cos 9 z -cf) bin. 

We assume that the neutrinos are collinear with the primary cosmic ray particles that 
produce them. The angle 9q u of the neutrino with respect to the shower axis, for approxi- 
mately half of the muon neutrinos, can be schematically written as: 9q u — 8^ ® 9^, where 
9 n is the angle between the parent meson (most of the time a charged pion) and the primary 
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particle, and 9 nu the angle between the meson and the neutrino; for electron neutrinos and 
for the other half of the muon neutrinos one has to consider a two-decay chain, and the 
angle between primary particle and neutrinos is 9 0u = 9^ © 9 nfJ , © 9^ u . Since the maximum 
p± kinematically allowed in a 7r ± ( y u ± ) decay is 30 (50) MeV, the dominant contribution 
to the neutrino angle comes from the transverse momentum of the parent meson (of order 
300 MeV): 



where we have used the fact that on average the neutrino has approximately one-quarter 
of the parent pion energy The angle 9q v is smaller that the angle 9g, v between the detected 
charged lepton and the neutrino, and its neglect does not introduce significant errors in the 
predictions of the geomagnetic effects on the angular distribution of the charged leptons. 

The next step is to treat the interaction of the neutrinos in the detector and find the 
direction of the produced leptons. For this purpose we use the quasi-elastic and single 
pion neutrino cross sections as calculated in Ref. |7| including corrections for nuclear target. 
For deep inelastic scattering we use the structure functions of Ref. ||. Neutral current 
interactions are neglected. The direction of each neutrino is chosen randomly within the 
bin and the direction of the outgoing electron or muon is then chosen randomly using the 
appropriate differential cross section. 

To provide realistic and relevant examples of the angular dependence, we consider two 
classes of events, applying cuts similar to those of Super-Kamiokande (SK) 0. As a low- 
energy sample, we use quasi-elastic simulated events in which electrons (muons) have mo- 
menta in the interval 0.1 < p e < 1.33 GeV/c (0.2 < p^ < 1.33 GeV/c). We compare this 
low energy sample to the 'single-ring' subset of the sub-GeV data of SK. As a high-energy 
sample, we compare all events with lepton momenta in the interval 1.33 < pi < 10 GeV 
with the multi-GeV data of SK. Our definitions are of course not precisely equivalent to the 
experimental classifications. For example, the single-ring events in the data include some 
multi-particle events in which only one is visible. Conversely, some quasi-elastic events 




(p ±j7r ) 300 MeV 



4.3° 



(4) 
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would be excluded from the single-ring sample because of an energetic recoil proton. We 
did check that both these contributions are small and that they do not significantly alter the 
angular distributions of leptons. Thus we believe that the simplified cuts we make on the 
Monte Carlo are adequate for our purposes to illustrate the expected angular dependence of 
the two categories of leptons. As a confirmation, we can compare the number of events in 
our cuts with the corresponding cuts in the SK data. We find 689 sub-GeV electrons and 
1050 sub-GeV muons (single ring only) as compared to 789 and 1185 in 25.5 kT-years that 
the Super-Kamiokande Collaboration report [f| from their simulation using the same neu- 
trino flux. For the same exposure, we find a total of 1034 multi-GeV leptons as compared 
to 1176 in the SK simulation. 

We show in Fig. la the distributions of neutrino energies that give rise to the two 
classes of events. The average neutrino energy for our sub-GeV muon sample os 0.8 GeV 
as compared to 5.7 GeV for the higher energy class. The order-of-magnitude difference 
in energy corresponds to a similar difference in L/E and makes the atmospheric neutrino 
beam a powerful probe of oscillations in an interval of parameter space with large mixing 
and 10~ 3 < Am 2 < 10~ 2 eV 2 . This follows from the well-known expression for survival of a 
neutrino flavor in vacuum, which in a two-neutrino example is 

Am 2 (eV 2 )L kr , " 



P vp -+v a = 1 - sin" 20 sin' 



1.27- 



(5) 



EGeV 

together with the large differences in pathlength between up-going and down-going neutri- 
nos. 

In Fig. |I|b we show the distribution of the angle between the detected muon and the 
parent neutrino. The average values of cos 9 for the sub- and multi-GeV samples are 0.53 
and 0.97 respectively. Normally only the charged lepton is detected, and because of the 
angle 6^ u with respect to the parent neutrino direction there is a smearing of the angular 
distribution of the neutrinos which is significant for the sub-GeV sample. Only detectors 
with high granularity can measure the recoiling nucleons (or more complicated hadronic 
final states) and reconstruct the neutrino energy and direction. Such a measurement is 
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potentially highly valuable in the search for neutrino oscillations. 

We illustrate the effect of the angular smearing in Fig. |2| by showing the azimuthal 
and zenith angle distributions for sub-GeV muons (solid lines) as compared to the same 
distributions for their parent neutrinos (dotted lines) at the location of SK. The azimuth is 
defined so that 0° corresponds to events from the north. Azimuthal angle increases counter- 
clockwise. The most prominent feature is the excess of events from the west. The cos8 z 
dependence for sub-GeV events at Kamioka is dominated by the high local geomagnetic 
cutoffs. We understand the slight excess of events from below (cos# z < 0) as arising from 
the fact that the local geomagnetic cutoffs, which affect the down-going events, are generally 
higher than the cutoffs averaged over the opposite hemisphere of the earth that regulate the 
up-going events. We note that both geomagnetic poles are below the horizon at Kamioka 
and therefore contribute to the relative excess of events for negative cos 9 Z . The depression of 
the neutrino flux near the horizon indicates that, for the sub-GeV events, the geomagnetic 
suppression more than compensates for the enhanced production of neutrinos from muon 
decay in this same angular region. 

In Fig. 3 we compare the azimuthal dependence for sub-GeV and multi-GeV events for 
four intervals of cos #2 of equal solid angle. The multi-GeV sample is sufficiently high in 
energy that the geomagnetic effects are much reduced. 



IV. RESULTS AND DISCUSSION 

Geomagnetic location is of great importance for the nature of the fluxes and angular 
distributions of low-energy events. We illustrate this in Fig. 4 by comparing the angular 
distributions expected in the absence of oscillations for sub-GeV muons at Kamioka with 



that expected at Soudan [JTO | or SNO ]TI|. The latter two experiments are near the north 
geomagnetic pole, so the flux of down-going events is significantly higher than at Kamioka. 
Moreover, the east-west effect is nearly absent for events coming from above. The sky-maps 
are in local coordinates with the local zenith at the top and the local nadir direction at the 
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bottom of each map. 

The angular-dependence of the neutrino events recently reported from Super- 
Kamiokande || has suggested several interpretations. Among these, the simplest possi- 



bilities are <-» v T and <-> v ste riie P^|~P^1- in Fig- 5 we show the expected zenith 



angle dependence for sub-GeV and multi-GeV neutrino induced muons at Kamioka (us- 
ing our definition). The solid line shows the result for no-oscillations. In the multi-GeV 
sample, the expected enhancement near the horizontal is clearly visible. The three broken 
lines show the results expected according to Eq. |5| for full mixing assuming «-> v T with 
Am 2 = 10 -2 , 10 -2 ' 5 , andlO -3 eV 2 . The distortion of the zenith angle distribution produced 
by neutrino oscillations depends on the oscillation parameters. Vertical up-going muons 
have the same suppression ~ 1 — - sin 2 29 because of the averaging of oscillations on the 
long pathlength. The shape of the suppression factor as a function of zenith angle depends 
strongly on Am 2 and is different in the two samples, reflecting the order of magnitude dif- 
ference in the typical energy of the neutrinos that give rise to the events. These features are 
potentially distinguishable with the future high statistics data of Super-Kamiokande. 

Whereas the shape of the zenith angle dependence strongly reflects assumptions about 
oscillations, the azimuthal dependence at fixed zenith is practically the same for all oscilla- 
tions scenarios. We show this in Fig. 6 for the sub-GeV muons (Kamioka, our definition). 
The four panels are for regions of equal solid angle of increasing zenith angle from the ver- 
tically down-going quadrant to the vertically up-going quadrant. Typical pathlengths in 
the two down-going quadrants are ~ 30 km and ~ 300 km, with large variations due to 
neutrino-lepton scattering angle as well as the relatively broad distributions of production 
height The up-going quadrants have pathlengths of order 10 4 km, with an admixture 
of shorter pathlengths near the horizontal direction. 

Table 1 gives a quantitative summary of the east-west effect for neutrino induced muons 
at Kamioka. The west/east ratio is ~ 1.35 for down-going and ~ 1.28 for up-going, sub- 
GeV muons. The ratio is ~ 1.10 for the multi-GeV muons. These ratios have a negligible 
dependence on the nature of the assumed oscillation. For this reason, study of the azimuthal 
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dependence of neutrino interactions should provide an important probe of the systematics of 
searches for neutrino oscillations with the atmospheric neutrino beam. In addition, study of 
the azimuthal and zenith angle dependence of electrons should be a sensitive test of whether 



(as suggested in Refs. |Il^ , [T4|j ) oscillations on terrestrial scales also involve electron neutrinos 
to some extent. 
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TABLES 

TABLE I. Average muon rates (in units (Kt yr sr) _1 ) in four solid angle quadrants. The rates 
are calculated in the absence of oscillations and for — > v T oscillations with maximal mixing and 
three values of Am 2 . 
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FIG. 1. a). Distributions of neutrino energies that give rise to the sub-GeV (solid line) and 
multi-GeV (dotted line) muon samples at Kamioka - see text for the definitions of the two groups 
of events in this calculation, b). Distribution of cos 6 (6 is the angle between the neutrino and 
muon direction) for the same two muon samples, same coding. 
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FIG. 2. a) Azimuthal distributions for the sub-GeV muons (see text) at Kamioka (solid line) 
and for their parent atmospheric neutrinos (dotted line) averaged over the zenith angle 6 Z . b) cos6 z 
distributions ( averaged over the azimuth angle (p ) for the same muon and neutrino samples. 
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FIG. 3. a) Azimuthal distributions for the sub-GeV muons (see text) at Kamioka for four bins 
in cos9 z : 1 to 0.5 (solid); 0.5 to 0. (dots); 0. to -0.5 (short dash); -0.5 to -1. (long dash), b) 
Azimuthal distributions for the multi-GeV muons at Kamioka, same coding. 
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FIG. 4. Two dimensional azimuthal and zenith angle distribution of the sub-GeV muons (see 
text) at Kamioka (top) and Soudan/SNO (bottom) in units of number of muons per kT.yr.sr. The 
top of the maps corresponds to the local zenith and the bottom to the local nadir. The north direction 
(<p = 0) corresponds to the edge of the map, south (ip = 180° ) to the vertical line in the middle 
with west (ip = 90°) and east (ip = 270°) to the left and right. Notice the smeared east-west effect 
in both maps. 
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FIG. 5. a) cos 6*2 distribution for the sub-GeV muons (see text) at Kamioka, averaged over 
azimuth - solid line. The other histograms are for — > v T oscillations with maximal mixing and 
Am 2 = 1CT 2 eV 2 (dots), l(T 2 - b eV 2 (dashes) and 1CT ?> eV 2 (dash-dash), b) The same distributions 
for multi-GeV muons, same coding. 
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FIG. 6. a) Azimuthal distributions for the sub-GeV muons (see text) at Kamioka for four bins 
in cos6 z : 1 to 0.5 (a); 0.5 to 0. (b); 0. to -0.5 (c); -0.5 to -1. (d) - solid lines. The other 
histograms show the distributions in the presence of — > v T oscillations for maximal mixing and 
and Am 2 = 10~ 2 eV 2 (dots), 1(T 2 - 5 eV 2 (dashes) and 10~ 3 eV 2 (dash-dash). 
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